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predictive  mode.  The  electron  density  model  relies  on  physical  Input 
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Section  1.0 


INTRODUCTION 


Most  quantitative  models  of  the  density  of  electrons  In  the  Ionosphere,  N^, 
constructed  to  date  have  not  been  usable  In  the  predictive  sense.  It  has  been 
realized  for  several  decades  that  the  variability  In  the  density  of  Ionospheric 
electrons  causes  problems  In  radio  communications.  The  nature  of  the  problem 
Is  thought  to  be  well  understood  physically.  Thus  the  availability  of  a quan- 
titative model  of  electron  density  that  can  be  used  to  help  predict  these 
variations  should  be  quite  useful.  The  Intent  of  this  contract  has  been  to 
develop  such  a model. 

From  the  beginning  It  was  realized  that  there  already  were  In  existence  several 
physical  models— that  Is,  models  which  attempt  to  describe  the  Ionosphere 
entirely  In  terms  of  first  principles.  There  are  at  present,  however,  several 
problems  with  such  models.  They  are  very  cumbersome  as  computer  codes.  None 
of  them  Is  global.  Rather  they  attempt  to  describe  a single  feature  of  the 
Ionosphere  or  region  In  latitude.  The  attitude  In  the  construction  of  the 
present  model  was  to  demand  from  the  outset  that  the  model  meet  several  require- 
ments. These  Included:  the  model  must  be  global— It  must  represent  N^  at  all 
latitudes  and  local  times.  It  must  be  very  fast  as  a computer  code  yet  accurately 
represent  all  gross  features  In  N^.  Most  Important,  the  model  must  be  constructed 
such  that  it  can  be  used  In  a predictive  mode.  Thus,  It  must  rely  on  physical 
Input  parameters  that  are  directly  measured  (the  Intensity  of  bands  of  the  solar 
electromagnetic  energy  spectrum  and  the  density  of  the  solar  wind  are  examples). 
Thisimplles  that  "secondary"  Indices  such  as  Kp  would  not  be  used  as  Inputs  to  the 
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model.  Such  Indices  are  averaged  over  long  periods  of  time  (a  few  hours)  and  are 
the  result  of  many  Interrelated  physical  events  In  the  near  earth  environment. 
There  are  other  "Inputs*  to  the  model  that  appear  In  the  form  of  constants  In  the 
equations  describing  the  Ionospheric  layers,  etc.  These  constants  must  be  deter- 
mined from  satellite  and  other  observational  data  sets  (the  "L"  value  for  the 
equatorial  anomaly  Is  an  example). 


A model  with  these  "primary"  (time  varying)  and  "secondary"  (fixed  program 
constants)  inputs  should  be  capable  of  predicting  when  the  proper  primary 
Inputs  are  used  In  real  time.  Currently  statistical  models  are  used  to  help 
with  radio  communication  problems  associated  with  variability  In  N^.  They  can 
be  used  only  to  suggest  seasonal  and  solar  cycle  trends  In  the  variability  In 
N and  are  therefore  of  little  use  In  predicting  the  day-to-day  and  hour-to- 
hour  changes  in  N^.  This  short  term  variability  In  N^,  however.  Is  of  great 
concern  for  comnunl cations.  While  the  present  model  does  not  address  such  fine 
structure  problems  as  "spread  F"  and  "sporadic  E"  It  should  be  useful  In  the 
study  and  prediction  of  the  variability  of  gross  spatial  structures  In  N^. 

The  model  features  In  Include  the  structure  of  the  E,  F^,  and  F2  layers,  the 
day-night  changes,  the  equatorial  and  seasonal  anomalies  and  dependences  on 
several  forms  of  solar  variability.  Thus,  solar  cycle  and  seasonal  variations 
find  their  way  Into  the  model  both  through  Its  primary  Input  parameters  and  In 
the  use  of  the  MDAC  model  of  atmospheric  density.  Both  electromagnetic  and 
corpuscular  energy  sources  are  used  to  provide  Information  for  model  Inputs. 

This  Inclusion  required  the  use  of  coordinate  transormatlons  between  solar- 
magnetos  pheric  and  geographic  and  geomagnetic.  Very  fast  analytical  transforma- 
tion codes  were  developed  for  this  purpose. 
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The  model  Is  complete  with  the  exception  of  “calibration."  Thus,  It  remains 
to  use  observational  data  as  Input  to  the  computer  code  and  then  compare 
calculated  values  of  with  those  obtained  from lonograms  at  the  appropriate 
times.  We  have  agreed  with  Naval  Research  Laboratory  personnel  that  the  month 
of  September  1977  will  be  used  as  the  epoch  for  the  primary  data  base.  It  was 
In  this  month  that  the  current  solar  cycle  began  In  earnest— there  are  periods 
of  both  quiet  and  moderate  solar  disturbance.  Hopefully,  several  data  sets 
from  S0LRA0>HI  will  be  used  as  Input.  The  necessary  lonograms  are  not  yet 
available  from  the  World  Data  Center,  but  we  are  told  there  will  be  sufficient 
data  within  a few  months  to  complete  the  model  calibration. 


A description  of  the  model  construction  and  computation  of  N^  based  on  the 


model  are  given  in  the  following  section. 
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Section  2.0 
MODEL  DESCRIPTION 


We  report  here  the  construction  of  the  model  In  terms  of  Its  primary  and  secondary 
constants  and  parameters.  The  model  specifically  Includes:  the  structure  of  the 
E,  , and  ?2  the  equatorial  anomaly,  the  winter  anomaly,  polar  cap 

structure,  semiannual  and  seasonal  dependences,  and  the  mid-latitude  trough. 

The  Input  parameters  to  the  model  are  chosen  such  that  satellite  data  can  be 
used  together  with  the  model  to  specify  Ng  In  real  time.  Hopefully,  with  data 
from  the  HELIOCENTRIC  satellite,  the  model  can  be  used  to  predict  Ionospheric 
features  controlled  by  magnetic  fields  and  solar  charged  particles.  In  part, 
the  predictive  capability  of  the  model  stems  from  Its  dependence  on  Inputs  from 
other  quantitative  environmental  models.  For  example,  the  location  of  the  mid- 
latitude trough  Is  largely  controlled  by  the  locations  (L  values)  of  the 
plasmapause  and  the  auroral  oval.  In  fact,  several  gross  Ionospheric  features 
are  controlled  by  the  geomagnetic  field.  The  MOAC  magnetic  field  model  provides 
the  necessary  magnetic  Inputs.  The  model  also  relies  heavily  on  the  MDAC  model  of 
upper  atmospheric  neutral  density,  N.  Several  Ionospheric  structures  and 
temporal  dependences  tie  directly  to  variations  In  atmospheric  density.  An 
example  Is  shown  In  Figure  1 where  the  semiannual  variation  In  at  the  F2  peak 
Is  shown.  This  variation  In  N^  Is  produced  entirely  by  the  semiannual  variation 
In  N as  described  by  the  MDAC  neutral  density  model.  The  structures  of  the  E and 
F-j  layers  are  also  determined  by  calculating  the  absorption  of  solar  electro- 
magnetic radiation  In  the  neutral  density  model.  A short  suimary  of  the  features 
of  the  MDAC  atmosoherlc  density  model  used  In  this  study  Is  given  below. 
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2.1  ATMOSPHERIC  DENSITY  MODEL 

The  MOAC  atmospheric  density  model  Includes  as  Inputs  both  the  solar  ultraviolet 
and  charged  particle  energy  sources.  The  model  Is  semi -empirical --It  Is  based  on 
available  satellite  data  and  our  present  understanding  of  both  the  UV  and  cor- 
puscular energy  sources.  It  offers  a global  description  of  the  atmosphere.  The 
model  takes  Into  account  the  dependence  of  atmospheric  density  on  the  value  of 
the  solar  flux  constant.  The  UV  and  corpuscular  effects  are  computed  In  separate 
coordinate  systms  and  their  contributions  added. 

Ever  since  an  accelerometer  was  first  flown  on  a polar  satellite  (Bruce,  1968), 

It  has  been  apparent  that  there  are  at  least  two  density  bulges  In  the  lower 
thermosphere  during  geomagnetically  quiet  times.  The  low-latitude  bulge  Is 
usually  attributed  to  heating  by  solar  UV  radiation  (although  other  energy 
sources  may  contribute  Importantly),  but  the  high-latitude  bulge  Is  produced  by 
particles  precipitating  Into  the  upper  atmosphere  through  the  dayside  cusps. 

The  solar  UV  heating  source  Is  best  described  In  geographic  coordinates.  Thus, 
Inputs  to  the  atmospheric  computer  subroutine  must  Include  the  universal  time 
and  the  time  of  year.  Local  time  Is  also  entered  as  a function  of  the  coordinates 
of  the  point  where  the  atmospheric  density  Is  to  be  determined.  The  dayside  cusp 
particles  are  constrained  to  precipitate  Into  the  atmosphere  along  magnetic  field 
lines.  The  cusp  Intersection  with  the  atmosphere  Is  nominally  15  degrees  below 
the  magnetic  dipole  axis  with  Its  longitude  center  on  the  magnetic  noon  meridian 
plane  as  defined  by  the  dipole  axis  and  the  sun-earth  line.  The  extent  of  this 
Intersection  In  longitude  Is  about  12  hours.  The  extent  of  the  cusp  In  latitude 
Is  several  degrees  (3-5).  The  region  of  the  atmosphere  actually  heated  by  these 
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particles  Is  of  course  much  larger  because  the  impact  energy  is  spread  out  by 
winds,  thermal  conduction,  and  possibly  by  gravity  waves.  The  corpuscular  con- 
tribution is  found  in  geomagnetic  coordinates  and  combined  with  the  UV  contribution. 

This  model  of  the  density  of  the  neutral  upper  atmosphere  is  functionally  very 
simple  and  computationally  very  fast.  Yet,  it  describes  most  of  the  known 
variations  in  the  atmosphere;  latitude,  longitude,  diurnal,  seasonal,  semi-annual, 
altitude,  solar  cycle  and  variations  in  particle  precipitation.  Thus  it  can  be 
used  to  predict  the  density  of  the  atmosphere  as  a function  of  time,  position, 
solar  and  magnetospheric  parameters.  Since  it  includes  the  heating  effect  of 
particles  precipitating  into  the  high  latitude  regions,  it  is  especially  valuable 
in  the  study  and  prediction  of  the  density  of  the  polar  atmosphere.  Furthermore, 
since  this  high  latitude  effect  is  parameterized,  it  can  be  used  to  predict 
increases  in  the  high  latitude  density  when  the  solar  wind  particle  flux 
increases  and  the  cusp  location  is  observed  to  change.  The  model  is  analytic 
and  thus  differentiable. 

The  three  dimensional  mercator  projection  in  Figure  2 shows  the  atmospheric 
density  at  400  km  during  summer  solstice  at  a universal  time  of  1600  hours.  The 
UV  density  peak  is  located  to  the  north  of  the  equator  and  at  about  1400  hours 

local  time.  The  particle  heating  peak  is  most  intense  on  the  subsolar  magnetic 

longitude  containing  the  magnetic  dipole.  Since  the  total  density  is  dependent 
on  the  particle  heating  effect  as  well  as  the  UV  heating,  the  northern  particle 

heating  peak  which  is  in  sunlight  is  much  more  pronounced  than  the  southern 

particle  heating  peak  which  is  far  into  darkness  during  the  summer  solstice 
night. 
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FIGURE  2.  Density  at  400  km  near  summer 
solstice  at  a universal  time  of 


2.2  LAYER  FORMATION 

In  addition  to  the  Inputs  from  the  MDAC  magnetic  field  and  atmospheric  neutral 
density,  the  model  of  N makes  use  of  our  knowledge  of  Ionospheric  layer  formation 
and  anomalies.  The  formulation  of  the  model  Is  discussed  below. 

Let  S be  the  energy  flux  at  location  Z and  S + dS  the  flux  at  location  l + di. 

Let  a be  the  absorption  cross  section  of  the  atoms  of  the  gas  and  N their  number 
density  (the  neutral  density  of  the  upper  atmosphere).  The  energy  absorbed,  dS, 

In  a cylinder  of  unit  cross  section  and  axis  parallel  to  the  direction  of  the 
Incident  beam  Is  given  by 
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If  n is  the  lon/electron  pair  production  efficiency  per  unit  energy,  then  q (the 
nuinber  if  ion/electron  pairs  produced  per  unit  volume)  is  given  by 

q{£)  ■ NotiS^  exp  (-x). 

The  above  equation  gives  the  production  of  electrons.  The  two  predominant  loss 
mechanisms  for  ion/electron  pairs  are  recombination  and  attachment.  For 
recombination  the  loss  equation  is 


-a 

e 


where  is  the  electron  density  the  positive  ion  density  and  a the 
recombination  coefficient. 


Thus,  at  equilibrium 


^ . q -c  . 0 


Then,  assuming  that  N 7^  N^, 


N<m  S exp  (-x)  -a  H ^ ■ 0 


Therefore,  N ■ [-— exp(-x)] 
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For  attachment  the  loss  equation  Is 


N^N 

In  discussing  attachment  processes  In  the  Ionosphere  It  Is  generally  assumed 
that  the  number  of  neutrals  Is  much  greater  than  the  Ions  so  that  when  Ionization 
occurs  the  neutral  density  does  not  change  significantly.  Thus 

Is  generally  used,  where  3 Is  the  attachment  coefficient. 

Then,  at  equilibrium 

dNe 

- q -6  Ng  » 0 and 

I 

Non  _ 

Ng  - — g—  exp  (-t) 

In  general  a and  3 may  vary  with  height  because  the  reactions  usually  Involve 
three  bodies  Instead  of  two. 
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2.2.1  The  E and  Layers 

It  has  been  observed  that  the  E and  the  layers  of  the  Ionosphere  can  be 
described  quite  well  using  the  recombination  equation.  Thus,  using  simple 
layer  theory  for  the  E and  F^  layers 


1 


Ng  - [N  a£  exp 

F, 

N - [N  Sp  Op  exp  (-Tp  )] 
e F^  F^  F, 


1/2 


E ^1  F 

where  and  are  the  electron  concentrations  for  the  E and  F^  layers,  S 
F] 

and  S are  the  solar  U.V.  flux  affecting  the  E and  F-j  layers,  and  Op  are 
the  recombination  coefficients,  N the  neutral  density,  and  Is  given  by 


fO 


N di  ^ Og 


N dz 


fO 


where  T = 


N dt  Is  the  total  atmospheric  cross  section  between  the  observation 


point  and  the  sun.  The  MDAC  atmospheric  density  model  (described  In  2.1)  Is  used 
In  the  evaluation  of  the  above  Integral. 


Then 


and 


Xe  - og  T 


rp  - Op  T 


where  and  Op^  are  the  absorption  coefficients  for  the  wavelengths  affecting 
the  E and  F^  layers. 
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These  equations  give  an  approximate  description  of  the  observed  structures  In 
the  E and  F-j  layers.  Changes  In  T depend  on  the  neutral  density  and  the  location 
of  the  sun.  Thus,  by  using  the  MDAC  neutral  density  model  (which  contains  most 
of  the  large  scale  variations  observed  In  the  neutral  atmosphere)  and  a simple 
sun-position  program,  the  temporal  variations  observed  In  the  E and  layers 
are  accurately  described. 


The  F^  layer  Is  essentially  turned  off  at  local  sunset  but  the  E layer,  while 
diminished  from  Its  daytime  strength,  persists  throughout  the  night.  Thus 

E I I 1/2 

Ng  (night)  = [N  exp  (-x^)] 

I 

where  Is  the  Intensity  of  scattered  light  and 

' 

“^E  “ ‘^E  N dt 
^0  00 

where  N dt  Is  the  optical  thickness  (Integrated  vertically  - not  toward  the 
sun).  00 

Profiles  of  Ng  at  noon  and  midnight  as  represented  by  the  model  are  shown  at 
the  equator  and  mid  latitudes  In  Figure  3. 

In  Figure  3 we  note  the  rather  large  unrealistic  dip  between  the  E and  F^ , and 
F^  and  Fg  layers.  This  large  decrease  In  electron  density  between  the  layers 
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Is  due  to  the  use  of  a single  wavelength  to  describe  each  layer.  This  single 
wavelength  description  Is,  of  course,  an  oversimplification.  The  electro- 
magnetic spectrum  from  the  sun  Is  quite  complex  and  several  wavelength  bands 
are  responsible  for  the  layer  formations. 

The  above  analysis  has  been  expended  to  Include  more  than  one  wavelength  In 
the  formation  of  the  layers.  The  only  time  consuming  calculation  Is  the 
Integral  which  calculates  the  optical  thickness.  This  calculated  optical 
thickness  can  be  used  for  all  wavelengths  having  different  absorption  and  Ionization 
coefficients  and  source  strengths  without  extensive  Increases  In  computer  time. 

Although  a large  number  of  wavelengths  can  easily  be  used.  It  was  found  that 
2 terms  for  E layer  and  3 terms  for  the  layer  provide  an  excellent  repre- 
sentation of  the  altitude  dependence.  The  model  currently  Incorporates  the 
following  wavelength  bands  In  the  parameter  list. 

E Layer 

Term  1 1026  “A  t 

2 977  “A  \ 

m 

10-170  “A 
170-796  “A 

796-911  ‘A  j 

i 

The  absorption  coefficients  associated  with  each  of  the  bands  has  been  adjusted 

s 

to  produce  the  electron  density  peak  at  the  correct  altitude  and  the  production  1 

coefficient  Is  adjusted  to  give  the  correct  magnitude  In  Intensity  at  the  peak.  I 


Layer 

Term  1 
2 
3 


i 
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atmospheric  density  model  (caused  by  the  slower  heating  effects)  Is  coupled 
with  the  rapid  response  to  changes  In  the  Ionizing  radiation. 

Figure  4 Is  a sample  altitude  profile  plot  which  combines  the  effects  of  the 
various  wavelengths  to  produce  a smoothly  varying  density  profile. 

2.2.2  The  F2  Layer 

The  structure  and  variability  In  the  F2  layer  Is  not  so  directly  controlled  by 
the  neutral  atmosphere.  Attempts  to  treat  the  F2  layer  using  U.V.  generated 
Ionization  and  electron  losses  which  are  a function  only  of  density  have  ended 
In  failure.  It  Is  generally  agreed  that  the  bottom  of  the  F2  layer  Is  defined 
as  the  region  where  the  electron  loss  rates  change  dramatically  (from  the  order 
of  seconds  to  the  order  of  hours).  Several  studies  were  undertaken  with  neutral 
constituents  and  electron  density  dependent  attachment  and  recombination  coeffi- 
cients. In  none  of  these  attempts  was  It  possible  to  directly  reproduce  the 
observed  seasonal,  semi-annual  or  solar  cycle  dependencies. 

An  empirical  function  was  therefore  developed  which  arbitrarily  (adjusted  to 
provide  a best  fit  at  some  point  In  time)  defines  a loss  term  varying  with 
height  and  an  Ionization  term  depending  only  on  the  ambient  density.  Further 
empirical  observations  suggested  that  the  normal  F2  layer  was  limited  to  the 

i 

i 
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ELECTRON  DENSITY  PROFILE 


Figure  4 


region  equatorward  of  the  latitude  of  the  plasmasphere.  Thus,  a function  was 


Introduced  which  limits  the  F2  layer  (except  for  the  winter  anomaly  peak)  to 
L values  ^ 4-6.  (The  L value  can  roughly  be  thought  of  as  the  equatorial 
extent  of  magnetic  field  lines  measured  from  the  center  of  the  earth  as  given 
In  earth  radii.  Thus  at  the  equator  L = 1.)  That  Is,  the  neutral  density 
dependent  F2  electron  values  are  multiplied  by  the  function. 


G(A)  » exp  (-Y  A"^) 


where  A » cos"^  \\l^  \ Is  the  magnetic  Invariant  latitude.  The  parameter  y 
was  adjusted  to  give  a best  fit  to  the  midnight  equatorward  edge  of  the  mid 
latitude  trough.  Thus,  the  parameter  y Is  a primary  parameter,  since  It  Is 
derived  from  the  size  of  the  plasmasphere.  To  properly  calibrate  the  parameter 
Y simultaneous  satellite  measurements  of  the  edge  of  the  plasmapause  and  measurements 
of  the  size  of  the  polar  cap  are  required.  (Note  however  that  the  model  will  pro- 
vide the  location  and  extent  of  the  trough  - that  measurement  Is  not  a model  Input.) 


The  F2  layer  Is  essentially  Independent  of  direct  solar  zenith  angle  control, 
and  depends  Instead  on  the  density  of  the  neutral  atmosphere  and  the  latitude 
extent  of  the  plasmasphere.  The  winter  anomaly, however,  does  appear  to  depend 
on  solar  zenith  angle.  It  appears  at  sunrise  and  rapidly  disappears  at  sunset. 
Much  of  the  average  structure  of  the  F2  region  as  given  by  the  model.  Including 
the  F2  peak.  Is  shown  In  Figures  5 and  6 In  terms  of  the  critical  frequency 
(which  relates  directly  to  N^).  Note  that  on  the  noon  meridian  some  of  the 
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structure  Is  produced  by  the  equatorial  anomaly  (around  the  12  Mhz  contour). 

The  magnetic  latitude  dependence  of  the  equatorial  anomaly  Is  Included  In  the 
model  and  Is  discussed  below. 

2.3  EQUATORIAL  ANOMALY 

The  equatorial  anomaly  is  thought  to  be  caused  by  an  upwelling  of  the  Ionosphere 
In  the  equatorial  region  produced  by  electric  fields.  It  may  be  organized  around 
a central  L value,  L^.  The  functional  form  of  the  anomaly  term  is 

d-exp  [-6(L-Ljj)^] 

where  the  amplitude,  d,  is  dependent  on  solar  cycle  as  well  as  local  time  and 
6,  which  defines  the  width  of  the  anomaly  peaks,  has  a local  time  and  solar 
cycle  dependence  and  a latitude  dependence.  At  a given  altitude  the  equatorial 
anomaly  is  strongest  along  the  magnetic  field  line  with  the  L value,  L^ 

(Lq  = 1.13).  Its  dependence  on  altitude  and  magnetic  latitude  as  given  by  the 
model  are  shown  In  Figure  7. 

2.4  WINTER  ANOMALY 

In  this  model  the  winter  anomaly  Is  treated  as  having  a completely  separate 
source  and  Its  contribution  to  the  electron  content  Is  assumed  to  be  additive 
to  all  other  electron  sources.  The  formalism  developed  for  the  layer  was 
also  used  to  describe  the  electrons  associated  with  the  winter  anomaly.  The 
winter  anomaly's  contribution  to  the  F2  layer,  N^  , during  the  winter  months 
near  solar  maximum  Is  therefore  written  as 
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MAGNETIC  LATITUDE 


A p 


F2u  1/2 

Ng  “ g(X,p)  [NSp  Op  exp(-Tp  )] 

® '^2W  '^2W 


where  g(X,p)  1s  a latitude  and  solar  cycle  dependent  function  which  determines 

the  extent  and  amplitude  of  the  winter  anomaly.  N Is  the  neutral  density,  Sp 

'’ZW 

Is  the  solar  flux  responsible  for  the  Ionization,  Op  Is  the  Ionization 

•■zw 

efficiency,  and  Xp  ■ Op  T,  where  Op  Is  the  absorption  coefficient  and  T 
•■zw  '’ZW  '’ZW 

the  optical  density  to  the  sun. 


In  general  there  Is  no  laye>*  over  the  polar  cap  In  tlie  above  formulation 
(except  for  the  winter  anomaly  which  at  times  extends  Into  the  polar  cap).  The 
layer,  however,  extends  to  higher  and  higher  altitudes  as  the  light-dark 
terminator  Is  reached.  | 

1 

2.5  POLAR  CAP  STRUCTURE  AND  MID-LATITUDE  TROUGH  ■ 

To  complete  the  polar  cap  structure  the  charged  particle  source  must  be  Included.  j 

These  high  latitude  sources  appear  most  often  at  auroral  latitudes  (at  L values 

greater  than  8.5-9).  At  night  the  particles  originate  from  the  tall  and  are  i 

1 

more  energetic  than  the  particles  which  continuously  precipitate  through  the 
dayside  cusp  regions. 

i 

The  configurations  of  the  auroral  E and  F layers  vary  extensively  In  amplitude  | 

1 

and  location.  For  sake  of  simplicity  a simple  Gaussian  In  latitude  has  been 

I 

used  centered  at  L • 10.5  with  a low  latitude  edge  of  8.5  at  midnight.  During 
daylight  hours  the  low  latitude  edge  moves  to  slightly  higher  latitudes  and  Is 
centered  on  the  dayside  cusp  field  line. 
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The  disappearance  of  the  F2  plasmasphere  associated  electron  density  (at  lati- 
tudes above  the  field  lines  passing  through  the  plasmapause)  together  with  the 
corpuscular  generated  auroral  electron  layer  forms  the  midlatitude  trough. 

During  daylight  hours,  especially  during  the  winter,  the  U.V.  produced  electron 
concentration  may  "fill  In"  the  trough.  The  variation  In  the  F2  peak  through 
the  trough  region  Is  shown  for  the  "average"  model  In  Figure  8.  Note  that  this 
Important  feature  In  Is  largely  under  magnetic  control.  That  Is,  as  the 
plasmasphere  expands  or  contracts,  the  low  latitude  edge  of  the  trough  moves 
respectively  to  higher  or  lower  latitudes.  Likewise  the  field  lines  along  which 
the  dayside  cusp  particles  precipitate  can  also  change  their  latitude  and  simul- 
taneously the  latitude  of  the  poleward  edge  of  the  trough.  The  comparison 
between  model  and  observations  of  the  motion  and  extent  of  the  midlatitude 
trough  should  be  one  of  the  first  quantitative  tests  made  of  the  model. 

2.6  SAMPLE  MODEL  RESULTS 

The  model  has  been  "calibrated"  using  average  relative  values  of  the  electro- 
magnetic flux  to  produce  average  profiles.  Some  of  these  profiles  for  the 
spring  equinox  are  shown  In  Figures  9 through  13.  Plots  9 to  11  show  the 
noon  time  density  profiles  for  various  latitudes  at  solar  minimum.  Plots  12 
and  13  are  examples  for  solar  maximum. 

In  order  to  better  represent  the  seasonal  and  solar  cycle  dependence,  plots  of 
the  critical  frequency  for  each  layer  versus  time  have  been  generated.  On 
these  plots  DAY  ■ 0 represents  January  1 at  solar  minimum.  The  electromagnetic 
flux  Is  then  allowed  to  change  In  a sinusoidal  manner  over  an  eleven  year  cycle. 
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Figure  9. 


ELECTRON  DENSITY  PROFILE 


Figure  10. 


089 


Figure 


Figures  14  to  18  show  this  solar  cycle  dependence  for  various  latitudes.  Note, 
for  example,  the  differences  between  northern  and  southern  latitudes.  They  are 
primarily  due  to  the  Interference  between  the  winter  anomaly  and  the  semiannual 
variation.  The  winter  anomaly  Is  hemisphere  dependent  while  the  semiannual 
variation  Is  not.  The  observed  noontime  critical  frequencies  at  Slough  are 
compared  with  model  calculations  In  Figure  19  and  shown  to  be  In  good  agreement. 
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Figure  17. 


Section  3.0 


MODEL  CALIBRATION 

The  description  of  the  model  in  the  previous  section  merely  provides  a glimpse 
at  its  "ground  state."  That  is,  the  figures  and  data  calculated  from  the  model 
are  based  on  "average"  input  parameters.  We  note  that  this  average  model  already 
contains  information  on  seasonal  variations,  dependence  on  solar  cycle,  and  the 


necessary  structure  in  latitude  and  altitude  in  order  to  accurately  describe 


"quiet"  features  in  N^.  We  believe  that  even  this  average  model  can  be  of  much 
value  since  it  is  global  and  has  an  associated  fast  computer  code. 


However,  the  purpose  for  the  construction  of  this  model  was  to  provide  a tool 
for  the  description  of  variability  in  N^.  This  is  only  possible  if  the  model 
possesses  the  proper  input  parameters  and  it  is  known  how  variations  in  each 
of  the  input  parameters  produce  change  in  N^.  The  model  input  parameters  have 
been  chosen  so  that  if  they  are  monitored  in  real  time  it  will  be  possible  to 
use  the  model  to  calculate  changes  in  N^.  They  include  information  on  various 
bands  in  the  solar  electromagnetic  energy  spectrum,  several  data  on  the  geo-^ 
magnetic  field,  and  a set  of  additional  variables  which  describe  other  solar 
and  Interplanetary  features. 


Note  that  magnetic  indicies  such  as  Kp  are  not  used.  There  are  two  reasons 
for  this.  First,  such  an  index  measures  a quantity  (the  range  in  a component 
of  the  earth's  surface  magnetic  field  over  a three  hour  interval  at  several 
polar  observatories)  that  is  controlled  by  many  geophysical  processes,  making 
it  difficult.  If  not  Impossible,  to  Isolate  causal  relationships.  Also,  these 
Indices  are  usually  not  available  on  a real  time  basis. 
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The  Input  parameters  used  in  the  model  then  are  chosen  so  that  they  can  be 
monitored  continuously  by  satellites  such  as  SOLRAD-HI  and  the  HELIOCENTRIC 
orbiter.  In  order  to  make  the  model  "responsive"  to  these  Inputs  It  Is  necessary 


1 

j 

1 

j 

1 


1 

j 

1 

to  calibrate  It— to  determine  the  response  of  the  single  model  output,  N^,  to 
changes  In  each  of  the  Input  parameters.  The  calibration,  of  necessity.  Involves 
two  types  of  data;  those  which  describe  variability  In  the  Input  parameters,  and 
proof  data— those  which  can  be  used  to  check  the  calculated  values  of  with 
the  observed  values  of  (associated  with  the  given  Input  parameter  data). 

The  principal  data  sets  to  be  used  for  this  calibration  are  from  the  SOLRAD-HI 
satellite  (to  describe  the  Input  parameters)  and  other  satellite  data  sets  as 
available.  The  proof  data  will  come  mostly  from  lonograms.  The  lonosonde 
stations  (latitude  and  longitude)  will  be  chosen  so  that  all  of  the  day-night, 

i 

latitudinal,  and  solar  activity  variations  In  the  model  can  be  calibrated.  In  j 

1 

agreement  with  the  NRL  people  and  NSF  (which  plays  the  role  of  US  coordinator  | 

i 

I 

for  the  International  Magnetospheric  Study)  we  will  concentrate  on  the  month  j 

of  September,  1977.  During  that  month  the  current  solar  cycle  began  and  there 
are  periods  of  both  quiet  and  moderately  disturbed  conditions.  Also,  the  IMS 
has  agreed  that  6 days  In  1977  will  be  analyzed  In  detail.  Three  of  them  are 
during  the  month  of  September.  We  therefore  expect  that  In  the  next  few  months 
there  will  be  a rather  large  amount  of  data  available  for  analysis.  The  lono- 
grams will  be  obtained  from  the  World  Data  Center  In  Boulder.  Their  people 
tell  us  that  our  requests  for  lonograms  can  be  honored  In  the  next  few  months 
as  they  receive  the  lonosonde  data  from  stations  of  Interest  to  us. 
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Once  calibrated,  the  model  will  hopefully  be  useful  both  for  specifying  gross 
features  In  and  for  the  actual  prediction  of  some  of  them.  Clearly  the 
Ionosphere  responds  quickly  to  changes  In  electromagnetic  radiation  from  the 
sun  and  It  Is  only  possible  to  predict  such  changes  by  first  predicting  changes 
In  solar  surface  features.  There  are  however,  several  aspects  of  Ionospheric 
behavior  that  are  Influenced  by  the  geomagnetic  field  and  the  precipitation  of 
charged  particles  Into  the  upper  atmosphere  at  auroral  and  polar  latitudes. 

Such  features  respond  to  changes  In  Interplanetary  parameters  such  as  the 
density  and  velocity  of  the  solar  wind  and  the  polarity  of  the  Interplanetary 
magnetic  field.  These  changes  can  be  measured  before  they  Influence  N^.  For 
such  changes  this  model,  when  calibrated,  should  constitute  a truely  predictive 


Section  4.0 


SUMMARY 


A global  model  of  the  density  of  Ionospheric  electrons,  N^,  has  been  developed. 

It  is  much  faster  as  a computer  code  than  any  of  the  existing  physical  models. 
Also,  it  is  not  statistical  but  rather  semi empiri cal . That  is,  a good  deal  of 
the  model  is  based  on  physical  understanding  of  ionospheric  features  but  other 
portions  are  based  on  observational  data.  The  model  has  been  constructed  such 
that  it  can  be  used  in  the  prediction  of  the  variability  of  N^.  It  is  intended 
that  the  model  inputs  are  all  based  on  data  observed  in  real  time  rather  than 
on  any  of  the  magnetic  indices.  The  model  can  of  course  be  used  In  other  than 
real  time  settings  and  with  minor  modification  used  with  magnetic  indices  as 
inputs.  We  expect  that  the  model  will  be  useful  to  the  Navy  in  its  continuing 
work  on  radio  communications  and  that  several  research  groups  will  find  it 
helpful  in  such  scientific  studies  as  the  understanding  and  description  of 
the  neutral  wind  problem  (which  takes  as  a starting  point  the  densities  of  both 
the  neutral  atmosphere  and  N^). 

The  calibration  of  the  model  with  satellite  and  ionosonde  data  should  be  complete 
in  the  next  several  months  and  the  model  available  for  use  by  the  end  of  the 
calendar  year.  We  would  be  happy  to  hear  any  advice  the  readers  of  this  report 
might  forward  to  us. 
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